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ABSTRACT 

ob 

, Several investigations in the study of cosmological structure formation use numerical simulations 

in both two and three dimensions. In this paper we address the subtle question of ambiguities 
in the nature of two dimensional gravity in an expanding background. We take a detailed and 
formal approach by deriving the equations describing gravity in {D + 1) dimensions using the 
action principle of Einstein. We then consider the Newtonian limit of these equations and finally 
obtain the necessary fluid equations required to describe structure formation. These equations 
are solved for the density perturbation in both the linearised form and in the spherical top 
hat model of nonlinear growth. We find that, when the special case of 13 = 2 is considered, 
^ . no structures can grow. We therefore conclude that, within the frame work of Einstein's 

CNJ I theory of gravity in (2 + 1) dimensions, formation of structures cannot take place. Finally, we 

indicate the different possible ways of getting around this difficulty so that growing structures 
ly-^ . can be obtained in two dimensional cosmological gravitational simulations and discuss their 

' implications. 
OO ' 

5 . 

. Subject headings: Cosmology : theory - dark matter, large scale structure of the Universe 
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Introduction 



The dominant paradigm for the generation of the observed large scale structure in the universe is based on 
the idea that the gravitational instability amplifies the initially small density perturbations. The equations 
' describing the growth of density perturbations in the highly nonlinear stage are analytically intractable and 

■ hence large scale numerical simulations arc resorted to for exploration of this regime. 

These N-Body simulations require large amount of computing resources (CPU, memory and storage space) 
if one is to get the requisite amount of dynamical range, i.e. good resolution in force and mass, large range 
in values of density etc. Time and resource constraints usually limit our ability to probe structure formation 
issues more deeply using computers, once the required resources are at the limits of technological feasibility. 
The key parameter which decides the feasibility level of numerical simulations is the size of a simulation, 
which — in turn — is characterised by: (i) The number of particles in the simulation volume, which is 
generally specified as , where D is the dimensionality of the simulation (usually 2 or 3) and (ii) The 
number of mesh points (M ) along any axis which determines the minimum length scale at which the results 
can be treated as reliable indicators of physical phenomena. In order to create a simulation volume that is a 
fair sample of the universe one needs about 10^ particles and in order to have a high enough force resolution 



one needs to increase the number of grid points adequately (for a review, see eg. Bagia and Padmanabhan 



1996| ). Let us suppose we have 160^ particles in three dimensions and our grid is 160 units on a side. 



Then, for the same amount of computational resources one can simulate a two dimensional situation with 
2048^ particles on a 2048^ grid (160'^ ~ 2048^). So, if we can extract useful (i.e. generalisable to the three 
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dimensional case) physical insights from results in two dimensions, then simulations of two dimensional 
gravity will be helpful. This hope has led to a large number of two dimensional simulations in the field of 
gravitational clustering (see eg Bagla et al. 199^ , Valinia et al. 1997 , Sathyaprakash et al. 1995 , Alimi et 



al. 1990| , ^handarin et al. 1989| ). 



There are three ways in which two dimensional gravity can be operationally defined and corresponding 
numerical simulations undertaken: (i) Consider a system of point particles in a three dimensional 
(expanding) background with the force of interaction being given by Newton's law of gravitation (i.e 
F (X 1/r^). The initial positions and velocities of the particles are such that they all lie in the same plane 
and all the velocities are in the plane i.e there are no velocity components orthogonal to the plane. This 
system will evolve with the particles being confined to the plane with clustering occurring in the plane. 
Thus, we have a two dimensional clustering scenario, (ii) Another system we can consider consists of 
infinite, thin 'needles' located parallel to each other. The mass elements in the 'needles' still interact 
through the 1/r^ force, but the interaction between 'needles' (obtained by summing over the mass elements) 
is given by a 1/r force. In this case as well, the background space expands uniformly in three dimensions. 
The two dimensional clustering that we study is the clustering of these 'needles', examined by taking a slice 
orthogonal to the 'needles', (iii) The third possibility involves writing down the Einstein's equations in two 
dimensions, finding the homogeneous and isotropic cosmological solution, taking the Newtonian limit (in 
which the potentials due to density perturbations and background metric can be superposed), finding the 
corresponding perfect fluid equations and solving them. In this case, we will have a background spacetime 
expanding in two dimensions unlike the other two cases. (There is yet another, fourth possibility, which can 
be defined only in an ad hoc manner. We will present this in the end.) 

The first case is not of much interest for cosmological simulations since the system is anisotropic, confined 
to a single plane and the clustering takes place in a specific plane only because the initial conditions were 
specifically selected to give this result. Hence we will not discuss it and it is mentioned here only for 
completeness. The way simulations in two dimensions are carried out usually is by simulating the second 
case and then defining the 'particles' as the intersection of the 'needles' with any plane orthogonal to them. 
In this case - as in the first case - the background spacetime expands in three dimensions (for a flat dust 
dominated universe the scale factor a{t) goes as t^/'^). The clustering that we observe and quantify in two 
dimensions, is basically the clustering of these needles in three dimensions. But this is also an anisotropic 
situation since the background spacetime is expanding in three dimensions. As an alternative we may try to 
write down the equations derived from Einstein's equations in two dimensions and examine how a system 
of particles interacting in a two dimensional expanding background spacetime is to be simulated. 

In the rest of this paper we shall examine the third alternative. We will take a very general approach by 
developing the formal theory of (_D + 1) gravity and considering Z? = 2 as a special case. (For some previous 



work on 2D gravity see eg. Deser, Jackiw and t'Hooft. 1984 and references cited therein 



The basic layout of the paper is as follows: In section (g) we first define the analogue of Einstein's gravity 
in (D + 1) dimensions, discuss the Newtonian limit and its corresponding Poisson's equation and then go 
on to analyse the Friedmann metric in (_D + 1) dimensions for a flat universe with dust. In section (H) we 
write down the D dimensional fluid equations and obtain the equation governing the density perturbations. 
This equation is then solved in the linear approximation and using the Spherical Top Hat model. Then in 
section (^) and section (|^) we specialise to the cases D — i and D = 2 respectively. Finally, in section (^ 
we summarise and discuss the implications of the results obtained in the earlier sections. 
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2. Formal {D + 1) dimensional gravity 

We start our analysis of {D + 1) dimensional (1 time dimension and D space dimensions) gravity from the 
action principle which we assume has the same form as that used in (3 + 1) dimensions. Using this action 
we construct the corresponding {D + 1) dimensional Einstein equations which will be subsequently used to 
study structure formation and spherical collapse. Thus, we begin with the action principle, 

5 = 5g = I d^^+^'x R^\ + / (1) 



2k{D) 

where iSg is the action for the gravitational field, is the action for the matter fields, g is the determinant 
of the metric tensor gik, R is the Ricci scalar, k{D) is a suitable constant which can be, in general, a 
function of D (when Z? = 3, k = SttG, G being the usual gravitational constant) and £m is the lagrangian 
density for the matter fields. The metric signature we adopt is (+, — ). We adopt the following 

convention regarding indices. Latin alphabets i,j,k... are used to represent (D + 1) dimensional indices 
which take on the values (0, 1, 2, . . . , Z?) while greek letters are used to denote D dimensional indices taking 
on the values (1,2,...,!)). Varying the total action S with respect to gik we obtain Einstein's equations. 



Gik = Rik ~ n 9ik R — —-r^Tik (2) 
2 



where Tik is the energy momentum tensor of the matter fields and is defined by 



1 rrirj. dW\g\C^) 




Gik is the Einstein tensor and Rik is the usual Ricci tensor. Note that the (1/2) that appears in Einstein's 
equations arises due to the square root in the term \/\g\ and has nothing to do with the dimension of the 
spacetime. 



We will use the above equations in the subsections to follow. In subsection (2.1), we will study the 
Newtonian limit of the metric tensor and then construct the corresponding Poisson equation that relates 
the Newtonian gravitational field to the matter density p. Then, in subsection ( |2.2| ), we analyse the 
Friedmann metric in (_D + 1) dimensions and the corresponding Newtonian limit of this metric is derived. 



2.1. Poisson equation in D dimensions 

In this section, we derive the Poisson equation relating the gravitational potential to the matter density 
p. We keep all factors of c since the Newtonian limit involves the limit c ^ oo. The analysis here follows 
closely the treatment in Landau and Lifshitz 1975 . Consider the metric. 

2 _ / 1 , 2(?!'\ 2 ,,2 ,,2 



ds^ = ( 1 + -T ) c^dr - dV (4) 



where is a function of space and time with dimensions of velocity square. The term is the D 
dimensional spatial line element given by the formula 



D 



dP^J2{dx"f (5) 
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We will now show that the metric written above is the Newtonian limit of Einstein's gravitational equations. 
We do this by showing that, in the Newtonian limit, the equation of motion of a particle follows Newton's 
force law with the force — mV<^... In relativistic mechanics, the motion of a particle of mass m is determined 
by the action function S 




ds — —mc cdti [l ^ — J- 5- (6) 



where v'^ is the square of the magnitude of the particle's velocity in D dimensions. In arriving at the second 
equality we have used the form of the metric in equation (^. In the limit c —^ 00, the action S can be 
approximated as 

„2 /"^W 1 I j ^ fdt (-mc^ + i?— 2 



S « -TOC^ dt[l+ ] ^ dt[ -mc^ + -mv^ - mcj) (7) 



The equation of motion for the particle can be immediately written down and we obtain 

m— = — mV(p (8) 

where v is the velocity vector in D space dimensions. Thus, Newton's force law is recovered in the 
non-relativistic limit and from this we conclude that the metric given in equation (^) is the Newtonian limit 
of Einstein's gravitational equations with (f> acting as the Newtonian gravitational potential. 

The relation between and the mass density p is found by taking the c — > 00 limit of Einstein's equations. 
This procedure, in (3 + 1) gravity, determines the constant k{D) since the Poisson equation is explicitly 
known. In other dimensions however, a definite criterion, like Gauss's law for example, must be imposed in 
order to determine k{D). We now consider the limit c — s- cx) of Einstein's equations in the following manner. 
First, we use the line element given in equation to calculate the Ricci tensor component i?oo 

where the summation convention has been invoked in the above equation and the sum over ji is only over 
the spatia/ dimensions. Then, using equation (^, we obtain. 



2k{D) 



^--tzttHtt^^ (10) 



where we have used the fact that gitg^^ = D + 1 and assumed D ^ \. Thus, Einstein's equations can be 
written in the equivalent form, 

R^k = ^ (r^k - ■T^-^9^kT] . (11) 



D-V 

where T is the trace of Tik. The energy momentum tensor of point particles is Tjfc = pc^UiUk where p is the 
mass density and Ui is the four velocity. Since, in the non-relativistic limit, the macroscopic motion is slow, 
the space components of Uj can be neglected and only the time component should be retained. Therefore, 
^0 = y/90Q and « for all fi. Consequently, only Tgo = 9oapc^ is non-zero. Substituting for Tik into 
equation ( pT| ) and using the expression in equation (|9|) for Rqo, we get, 

1^..,,, fD-2\{ 2(f>\ k(D) 2 fD--2\K(D) , , 

That is, 



D~2 
D ~ 1 



n(D)p (13) 
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where is the usual Laplacian operator in D dimensions. This equation is the Poisson equation in D 
dimensions. Note that when substituting for the value of Rqq from equation (|^), we neglected the first term 
in comparison with the second since the former is of order while the latter is only of order . 



2.2. Priedmann Universe in (1 + D) dimensions 

Let us next consider the maximally symmetric Robertson- Walker metric in (_D + 1) dimensions, specialising 
to flat space with A: = (we set c = 1 in this and in subsequent sections) , 

ds^ = dt^ -a^it)dP (14) 

where a{t) is the scale factor and dP is the D dimensional line element given in equation (H). Calculating 
the components of the Einstein tensor, we obtain, 

D{D -l)d^ 

LfOO — 



2a2 

Gn=G22 = ■■■ = Gdd = {l-D)ad+(^l-j^{D-l)d^ (15) 

where d stands for da(t)/dt and similarly d is the second derivative of a{t) with respect to time. All 
the other components are zero. For consistency, the energy momentum tensor must have the form 
r\. — diag(/9, —p, ~p, -—p, . . .) where p is the matter density and p is the pressure. 



Substituting in Einstein's equations, we obtain, 

D{D -l)d^ 



2a? 

d D - 2 a2 k{D)p 



n{D)p (16) 

(17) 



a 2 a? D-l 

The above two equations, together with the equation of state in the form p = p{p) completely specify 
the system. Solving these three equations, we can determine a(t), p(t) and subsequently p{t). Combining 
equations (16p7|), we get the single equation. 

We now specialise to the case of pressureless dust with the equation of state p = 0. Using the principle of 
conservation of energy and momentum expressed by the relation 

T\k = (19) 

we derive the following relation, 

, , — I V li/K i] - 
\9\ 

Noting that the only non-zero component of T\ is Tq = p wc finally get 

pa^ — constant — Ci (21) 



rr.k 1 d f /]— r„fc\ 1 dgkl r^k 
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Substituting the above relation into equation (|T^, and solving for a(t) and subsequently for p(t), we obtain 
the solutions, 

Let us next consider the Newtonian limit of the Friedmann metric. This limit is important because the 
length scales of interest in structure formation are small compared to the Hubble radius and the velocities 
in the system are also much smaller than c. This permits us to study the formation of large scale structures 
in the universe in a Newtonian framework where the effective potential due to the expanding background 
universe, $frw, and the potential due to the density perturbations, (p, can be simply superposed. In order 
to obtain <&frW: we first recast the Friedmann metric in equation (^4|) into the more convenient form 

ds'^ = dt^ - a^{t) {dX^ + X^dn^) (23) 

where X is the radial distance in D dimensions and is the corresponding solid angle. We then apply the 
transformations (see Padmanabhan 1996, pg 80,346) 

r = Xa{t) , T = t - to + ^adX^ + OiX"^) (24) 

where only terms up to quadratic in X are retained. Direct calculations, correct upto this order, transforms 
the Friedmann line element to the form, 

ds^ « (^1 - dT^ - dr^ - r^dn^ (25) 

which upon comparison with the metric in equation (^) gives the equivalent Newtonian potential $frw in 
D dimensions as 

*FRW ^-l-r^ (26) 
2 a 

We will now use the results developed in the last two subsections to study structure formation and spherical 
collapse using the STH model. 



3. Structure formation in D dimensions 

Having determined the form of the Poisson equation in the Newtonian limit and analysed the Friedmann 
equations in (D + 1) dimensions, we proceed to derive the equation for the growth of inhomogeneities in the 
expanding universe. After this we consider a specific model, the STH model, to study spherical collapse of 
matter. 



3.1. Equation for density perturbations in D dimensions 

Let us assume that matter in the universe is a perfect, pressureless fluid with density pm and flow velocity 
U. We can formally write down the D dimensional fluid equations describing a perfect fluid in an external 
potential field <i>tot in a proper coordinate system labelled by the D dimensional vector r. Therefore, we 
have, 



d Pa 

dt 



V7r-(/5mU) = (27) 
(U-Vr)U = -Vr$tot (28) 
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where equation (|2^) is the usual continuity equation while equation ( |28D is the Euler equation for the fluid. 
The potential in equation (p8|), $tot, is the total external Newtonian potential 

«'tot = «'FRW + V (29) 

where $frw is the background potential associated with the smooth background matter density pbm and is 
given in equation (p6|) while ip is the potential caused by density perturbations (pm — Pbm)- The potential 



Lp satisfies the Poisson equation given in equation (13). Thus, 



V?(p = «(^)(Pm - Pbm) = (§7:1) <D)p^^5 (30) 

where 5 is the density contrast defined by 

Pbm 

We now transform to comoving coordinates defined by x = r/a{t) and define the peculiar velocity v by the 
relation 

U = i/(t)r + V = ax + V (32) 
where v = ax and H{t) = {a/a). Then, equation ( p7| ) and equation ( p8| ) become 

^^""^ ' Z?i7p,„ + iv„ • (p,„v) = (33) 



at a 

dv\ 1 1 

— + Fv+-(v- Vx)v= — Vx(^ (34) 
ot J ^ a a 

where we have used equation ( p6| ) to substitute for $frw- Similarly, in co-moving co-ordinates, equation ( |30| ) 
reduces to 



Vi^ = ( ] <D)a'pkmS (35) 

Using pni = Pbm(l + S), transforming the time variable from t to a(t) and defining a new velocity variable u 

by 

dx V , , 

u = — - — 36 

da aa 

we can obtain equations for S{a) and u(a). Therefore, using equation ( |2l| ) and performing the 
transformations, equation (^) and equation (|3j) further reduce to, 

|^+Vx-[u(l + <5)]-0 (37) 
oa 

a + 2— u + a2(u- Vx)u = -W^p (38 

oa \ a J 

Now, we use the Friedmann equations in equations ( p^ , p7| ) with p replaced by pbm and with p = to 
substitute for d in the above equation. Further, we define a new potential \l/ by the relation 

V 6 - D / K[D)pi,^a-^ 
so that, upon using equation (|35|), one obtains. 
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where all reference to k{D) has disappeared. Hence the final system of equations we need to tackle are, 

1^ + Vx • [u(l + <5)] = (41) 
^" + (u • Vx)u = -^—^A [Vx* + u] (42) 



da 2a 



where A is given by the relation 



f 2k{D) \ a2 pbm(t) _ f D{D-l) \ 

For the fc = universe, we will set A = 1. 

To proceed further and determine the equation satisfied by 5, we decompose the term daU(i (where up is 
the /3th covariant component of the vector u and da is short for d/dx") as 

daUi3^(Jap + ^ap + -^5apO a,/3= (1,2,...,D) (44) 

where (Tq,/3 is the traceless, symmetric shear tensor, ^lap is the antisymmetric rotation tensor, 9 is the (trace) 
expansion and 5ap is the Kronecker delta symbol. Then, equation ( ^l]) and equation (^) are combined 
by taking the divergence of equation ( p2[ ) and using the above decomposition of daUp to obtain a single 
equation for 5. Straight forward algebra gives. 



(P5 (Q-D\d5 (D{D-2)\ fD + l\ 1 ( dS^ 



d^?n-^j^-l-^J'^' + '^^l^J(TT^UJ +(l + ^)(a^-2f^^) (45) 

where cr^ = aapcr"^ and fi^ = {l/2)Qap^"^ ■ This equation is the full non-linear equation for S. Apart 
from the obvious nonlinear terms containing (5^ and {dS/da)'^, the term (1 + 5){a^ — 2r2^), which is the 
contribution from the shear and rotation, is also non-linear. The non-linear terms in S in the above equation 
render the equation unsolvable in general. Ignoring these non-linear terms to a first approximation, we can 
get a linear equation for (5(a), 



d^S /6-D\ dS (D{D-2 



(5 = 0. (46) 



da^ ' \ 2a J da \ 2a 
Assuming a power law solution for delta in the form S (x a^, we get, 

p = ± \\/%D^ - 2AD + 16 (47) 

4 4 ^ ' 

as the required values for p. Notice that 5 has a growing mode as well as a decaying mode in general. The 
above solutions hold for all values of 13 > 1 in the linear regime. 

Though the full non-linear equation is not solvable, by neglecting the contribution from the shear and 
rotation terms and by using a suitable ansatz for 5, the resulting non-linear equation can be solved. We 
proceed to do this within the framework of the STH model in the next section. 



3.2. The Spherical Top Hat (STH) model 



In the STH (spherical collapse) model, we assume spherical symmetry by neglecting the shear and rotation 
terms in the equation for 6. With this assumption the 5 equation can be exactly solved. 
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Transforming equation ( jiq ) by changing the independent variable back to t, dropping the rotation and 
shear terms and using the Friedniann equations given in equations ( |l^ , |l7| ), we get, 

d^S ddS fD + l\ 1 fdSV fD-2\ ^. ^. 

We now define a function R{t) by the relation 

l + 5= — = (49) 

where Cd ~ 27r^/^/(_Dr[_D/2]) is the volume of a unit sphere in D dimensions introduced for later 
convenience and M is a constant. The expression for S above can be rewritten using the relation 
Pbma^ = Poo-o from equation (pl|): 



l + S = 



M 



D 



= ^J^ (50) 



where po and oq are the matter density and scale factor at some (arbitrarily chosen) "present" epoch tQ. 
Substituting equation (|o|) in equation (^), we get an equation for the growth of R{t) as, 

d'^R _ D~2 k{D) M 
^ ^ D{D - I) Cd 

[As an aside we may note that if the universe contains matter or fields with equations of state other than 
p — 0, the equation for R{t) becomes 

d'^R D-2 k(D) M k(D) , , 

-di^ = -DiD^)-^R^ - I^((^-^)^ + ^^)-^ (^2) 
where the term {{D — 2)p + Dp)rest comes from the smoothly distributed component with p ^ 0.] 

^From the form of the equation of motion of R{t) we can give the following interpretation. Since the entire 
system considered above is spherically symmetric, we interpret R as the radius of a D dimensional spherical 
region containing a mass M. The equation of motion of R determines the motion of the surface of this 
region. In general, a spherical overdense region will be expected to initially expand because of the expansion 
of the background universe till the excess gravitational force due to the overdensity of enclosed matter stops 
the expansion and causes the region to collapse back on itself. We will discuss the cases D = 3 and D = 2 
in the subsequent sections and determine the differences in the behaviour of the growth of inhomogeneities. 



4. Summary of standard results in 3-dimensions 

When D = 3, all the standard equations are recovered. First the Poisson equation satisfied by the 
Newtonian gravitational potential given by equation (p^) reduces to the standard form, 

= AnGp (53) 

where we have defined G by relating it to k(3) by k(3) = 8nG. Similarly, equations (^ij) and equation ( ^^ 
reduce to (with A= 1), 

^ + V^-[u{l + d)]=0 (54) 
du , ^ , 3 ^ 

— + u- Vx u= -— Vx* + u 55 

oa 2a 
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while the equation for ^' becomes, 

V^* - -. (56) 
a 

In a similar manner, the S equation reduces to, 

(PS 3 dS 3 ^, 4 /d(5^ ^ 



da^ 2a da 2a^ ^ ' 3(1 + 5) \da ) ^ > ^ > 

and the solutions to the linear perturbation equation which is obtained by dropping the nonlinear terms 
and the {a^ — term, are 

(Scxaf, P=l,-^ (58) 

which are well known. The STH model for Z? = 3 also reduces to the standard form 

d^R GM AnG , , 

— (/' + 3p),,,i? (59) 



which is again well known (see Padmanabhan 1996). Therefore, it is seen that the full D dimensional 
equations reduce to the correct equations in three dimensions. Now we will go on to discuss the important 
case of Z? = 2. 



5. Two dimensional Gravity 

If we naively consider the limit D — > 2 in the D dimensional equations, assuming that k(-D) is finite in this 
limit, we obtain the following results. First, the Poisson equation ( p^ reduces to 

V^^ = (60) 

The above result shows that in two dimensions, the gravitational potential does not couple to the matter 
density p. In structure formation, this means that inhomogeneities cannot grow since the perturbed 
potential Lp is not related to 8 at all. The second interesting result is that the background Newtonian 
potential ^frw vanishes. This occurs because, referring back to equation (p^), a = for pressureless dust 
and hence the background potential is zero. Further, the b equation reduces to 



d^b 2dS 3 fdS'^^ 



^) +{1 + S){a'^ -2n'). (61) 



da^ a da 2{1 + 5) \da 
Linearising the equation as before by dropping the (cr^ — 2il'^) and {dS/da)'^ terms we obtain 

d'^S 2dS „ 

y^ + -:r = o. (62) 

da a da 

The solutions to the linearised equation are 

d(xaP, p = 0,-l. (63) 

Thus only a constant or the decaying mode is present. This is consistent with the result that the perturbed 
gravitational potential does not couple to S. If one considers the STH model, it is easy to see that the 
growth equation for R{t) reduces to 

d^R 

— = -K(2)p„sti? (64) 
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For, prcst — 0, the solution to the above equation is just R(t) = B\t + where _Bi, B2 arc constants. This 
is to be expected since there is no gravitational force which can lead to clustering and as a consequence the 
radius simply grows with time just like the background universe. Thus, if is finite in the limit ^ 2, 
it is not possible to have gravitational clustering that can grow with time. 

We can, however, try some alternative approaches to examine whether it is possible to have a consistent 
physical picture of growing structures for two dimensional gravity. 

One possibility is that instead of assuming k(-D) to be finite, let us assume that the expression (k(_D)(Z? — 2)) 
remains finite when D — > 2. This finite value can be fixed, for example, by invoking Gauss's theorem in D 
dimensions. This gives 

Thus, k{D) — + 00 when D 2, but the Poisson equation acquires the form 

= 27rGp (66) 

This is, of course, the same form which is obtained by applying Gauss's law in two dimensions. Hence, as in 
three dimensions, the gravitational potential is determined by the matter density and thus inhomogeneities 
can in principle grow. There are, however, difficulties with this approach. To begin with, the constant 
factor / {2k{D)) in the action S in cquation(|^) vanishes for D = 2. But this is not too serious a problem. 
The gravitational part of the action certainly vanishes but because only the variations about the action 
are of significance this difficulty can be ignored. But a more serious problem arises when the solutions 
to the Friedmann equation are considered. The solutions for a{t) and p{t) in D dimensions are given in 
equation (p2|). Using equation (|65|), these reduce to, 

'm'i ."^":^S °f'" ■ .w = fto-=<£^m<- (67) 



{D-2)T[D/2]J ' J Dtt^/^G 

When Z) — > 2 then, a ^ 00 and p — > irrespective of the dependence on t. This implies that one cannot 
solve the equations describing the growth of structure in a consistent and non-singular way. Hence, 
we conclude that it is not possible to have a theoretical formulation of two dimensional gravity as the 
Newtonian limit to Einstein's equations in two dimensions. 

An alternative that remains is to use the Newtonian fluid equations in D dimensions directly and rewrite 
them for an expanding background with an arbitrary scale factor a{t). Note that a{t) is not obtained from 
the Friedmann equations and is completely arbitrary. We can superpose the potentials for the background 
universe and the perturbations in this case as before. The further assumptions we need to make are (i) the 
potential of the background universe $bg is of the form 

and (ii) the Poisson equation is given by 

V^^ = n{D)p^^5 (69) 

where k{D) = 2tt^ /'^G /T[D /2]. This form of k{D) is obtained from the use of Gauss's law in D dimensions. 
We also need to specify how the background density pbm depends on time. In analogy with the usual 
Friedmann equations, we will assume phmO-^ = Ci where Ci is a constant. This gives an equation for S 
with an arbitrary scale factor a(t) as 

->^iD)C,-j^&{l + 5) =[—-]-—[-] +{l + 5){a^-2if). (70) 



da^ \ ab? I da a^a? \ D J 1 + d \da J 
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The above equation can be solved in any dimension D if the form of a{t) is given. This gives us a 
non-singular way to analyse growth of structures in D dimensions, including the case D — 2. But in three 
dimensions we observe that the above equation does not correctly reduce to equation (|5^). We may obtain 
the correct equation in three dimensions by making an additional ansatz, namely, that 



= Pbm- (71) 



With this, equation (It^) reduces to 



Notice that the above equation differs from the earlier equation for S in D dimensions, equation (p5|), in 
that there is a factor of {D — 2) missing from the coefficient of the S{1 + S) term. Therefore, the above 



equation does correctly reduce to equation (57) since {D — 2) equals unity when D = 3. When D = 2, 
equation (^ gives 

-S{l + d)^— K 3- + {1 + S){a^ - 2n^) (73) 



da^ a da a^ 2(1 + (5) \da 

On linearising this equation by dropping the (ct^ — 2il^) term and the nonlinear terms and solving it wc get 
both a growing mode as well as a decaying mode for 5. The solutions are 

5cxa«, g=(-l±\/5)/2 (74) 

(It is interesting to note that one of the power law exponents is the golden ratio). 
The Spherical Top Hat (STH) equation in this case turns out to be 

f, GM IR 

^ = -^ + 2^^ 

where M is the constant mass inside a 'spherical' shell of radius R. This equation, unfortunately, has no 
simple analytic solution. 

While this procedure leads to nontrivial results, it has many ad hoc assumptions and cannot be obtained by 
taking appropriate limits of Einstein's theory in a systematic manner. Consequently it cannot be applied 
to numerical investigations of 2 dimensional gravity with the confidence that the results will have some 
implications for the three dimensional case. 



6. Conclusions 

In this paper we have analysed the case of two dimensional gravitational clustering starting from a 
formulation of the D dimensional Einstein's equations and taking the proper limits. The system of 
equations thus arrived at for a (D + 1) dimensional universe has been shown to reduce to the correct 
equations in three dimensions. But when the I? — > 2 limit of these equations is taken, we are forced to 
conclude that irrespective of the value of k{D), a consistent two dimensional gravity theory in a cosmological 
context that supports growth of structures cannot be constructed. 
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If k(2) is assumed to be finite, we observe that the coefficient in Poisson's equation goes to zero thus 
decouphng the potential from the density. This implies that perturbations do not grow but decay in 
time due to the expansion of the background spacetime. The alternative which is obtained by using the 
expression for k{D) given by equation (^) gives rise to solutions for the scale factor which are singular and 
therefore unacceptable. 

We have discussed all the ways in which two dimensional gravity may be simulated including an ad hoc 
procedure without a strong foundation which can give non-singular results as far as structure formation 
scenarios in two dimensions are concerned. The results presented in the paper leads us to conclude that 
the only way to do a numerical simulation of two dimensional gravity is to simulate infinite 'needles' in 
a background spacetime expanding in three dimensions and consider the 'particles' in the system to be 
intersections of the 'needles' with any plane orthogonal to them. 
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